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Abstract 

We consider the transverse-momentum distribution of heavy flavours in 
photon-hadron colhsions. We present a formahsm in which large transverse- 
momentum logarithms are resummed to the next-to-leading level, and mass 
effects are included exactly up to order aemOs) so as to retain predictivity 
at both small and large transverse momenta. Phenomenological applications 
relevant to charm photoproduction at HERA are given. 
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1 Introduction 



This work deals with the computation of the transverse momentum distribution in 
heavy flavour photoproduction. At present, fixed order (FO) calculations are avail- 
able, including NLO (Next-to-Leading-Order) corrections fl], ^. Furthermore, for 
very large transverse momenta, the so-called fragmentation function (or resummed) 
formalism, that allows to resum enhanced terms of order aemtts(tts logPi/'^^)* (which 
we call leading-logarithmic terms, or LL), plus terms of order aem^s (o^s logpT/^^)* 
(next-to- leading logarithmic terms, or NLL) has been developed in ref. 0], building 
upon ref. [^]. This approach has however the drawback that it is essentially a "mass- 
less" formalism, in the sense that it does not include contributions to the cross section 
that are suppressed by powers of m/pT- 

Several HI and ZEUS [^] results are presented in comparison either with calcu- 
lations performed in the resummed (massless) approach (RS), or with the fixed-order 
NLO calculation. It is thus important to provide a framework for computing a heavy- 
flavour cross section which is accurate in both the large and the small transverse 
momentum regions. A method for performing the merging of the FO 0, ^ ^ |lOl 
and RS [TTI calculations in the hadroproduction case was developed in ref. [jT^. In 
the hadroproduction context it was found that the mass corrections are positive and 
large, a result somewhat contrary to the intuitive belief that masses reduce the phase 
space, and thus the cross sections. 



The aim of the present work is to extend the formalism of ref. [T^ to the pho- 
toproduction case. This extension is not a straightforward one. In fact, it is well 
known that a generic photoproduction cross section has to be written as the sum of 
two components {pointlike and hadronic, also called direct and resolved respectively). 
For the transverse momentum spectrum of a heavy quark we write 
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where pnt stands for pointlike and hdr for hadronic; y is the rapidity of the heavy 
quark in the laboratory frame, and 
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The sums run over parton flavours, and are the four-momenta of the incoming 
photon and hadron respectively, and rfa^j, d&ij are the subtracted partonic cross 
sections. The pointhke and hadronic components of eqs. ( [1.2|) and (|1.3|) are strictly 
related beyond the leading order in perturbation theory: none of them is a physical 
quantity, only their sum (eq. ( |1 . 1| ) ) is measurable. This is the origin of the most 



serious problem we face when extending the formalism of ref. [0 to the present case. 
In this work we shall follow the strategy adopted in that paper, but we shall point 
out the major differences with respect to it, due to the problems inherent to eq. ( |1 . 1| ) . 

In perturbation theory at next-to- leading order (i.e., the highest accuracy reached 
so far in the computation of heavy flavour cross sections), we have the expansions 

da.yj = aemC(sda!^j + aemOtlda''^^ , (1.4) 
dcjij = a^da^^^ + a^da^f . (1.5) 

However, the parton densities in the photon Pj^'^^ behave asymptotically (i.e., at large 
scales) as agm/ag. Thus, at least at the formal level, the perturbative expansions 
of the pointlike and hadronic components of eqs. ( [1.2|) and ( |1.3|) are both series in 



aem«s, as in eq. (pT^). This allows us to simplify substantially our presentation; in 
what follows, we shall write the physical cross section of eq. ( p..l|) as an expansion 
in aemCtg- The reader must keep in mind that, in doing this, we are not referring 
to the pointlike component only, but to the observable that is actually measured 
in experiments. When we shall deal either with the pointlike or with the hadronic 
component only, we shall indicate it explicitly. 

Having clarified this point, we proceed in our program of extending the formalism 



of ref. |12| to photoproduction reactions. This means that we shall implement a 



computation with the following features: 

• All terms of order aem«s and aem«s included exactly, including mass effects; 

• All terms of order aem«s (^s IosPt/^^)* and aemttg (o^s logPT/"^)* are included, 
with the possible exception of terms that are suppressed by powers of m/prj.. 

To be more specific, let us write schematically the result of the NLO calculation of 
the photoproduction cross section as 

, 9 = A{m)aemOis + B{m)aemal + C(aemas) . (1.6) 
dydp^ 
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The explicit dependence of A and B upon (the centre-of-mass energy), y, p^, and 
the factorization/renormahzation scale /i is not indicated, and = as(/^)- The NLL 
resummed cross section is given by 



oo oo 



= aemas"^ ai{a Jog n/my + aemaj'^biia Jog n/m) 



dydp, .=0 .=0 

+0{a,^al{a, log/i/m)^ + 0{a,^a, x PST) , (1.7) 

where PST stands for terms suppressed by powers of mj'p^ in the large-px limit 
(possibly with further powers of mass logarithms). The coefficients and hi depend 
upon Ecmi Vi Pt and /i. If /i ~ p^, they do not contain large logarithms of the order of 
logpT^/m. The only large logarithms are the ones explicitly indicated. Our approach 
combines the results of eqs. ( |1.6| ) and ( |1.7| ), giving 

da 



dy dp] 

f oo oo \ 

'^ai^ajogiijmy ^ a^^al'^Y^hiiaJogiijmf j x G{m,p^) 

=2 i=\ ) 

^0{a,^al{aJog^lmy) + 0(aem«s x PST) , (1.8) 

where the function G{m,prj,) is quite arbitrary, except that it must be a smooth 
function (also in the p^^O limit), and that it must approach one when m/px— 
up to terms suppressed by powers of m/pT- Observe that the sums now start from 
i = 2 and i = I, respectively, in order to avoid double counting. Thus, this formalism 
contains all the information coming from the fixed-order NLO calculation and from 
the NLL resummed calculation. The arbitrariness in the function G arises from the 
fact that we do not know the structure of power-suppressed terms in the higher orders 
of the NLL resummed calculation. The choice of the function G only affects terms of 
order aemC^sy so far unknown. 

Before turning to the practical implementation of eq. ( |1.8[ ), we stress that it is 
important that both the RS and the FO approaches are expressed in the same renor- 
malization scheme. The commonly used FO approach uses a renormalization and 
factorization scheme in which the heavy flavour is treated as heavy. Thus, if we are 
dealing with charm, we use as of 3 light flavours as our running coupling constant, 
and the appropriate structure functions should not include the charm quark in the 
evolution. The RS approach, on the other hand, also includes the heavy flavour as 
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an active, light degree of freedom. This problem can be easily overcome by a sim- 
ple change of scheme in the FO calculation. Section |^ contains the details of this 
procedure. 

Once this is done, the FO calculation matches exactly the terms up to order 
ttemtts resummed approach, in the limit where power-suppressed mass terms 

are negligible. In order to subtract from the RS result the terms already present in 
the FO, we must provide an approximation to the latter where terms suppressed by 
powers of the mass are dropped. We shall call FOMO this "massless limit". In the 
simplified notation of eqs. (|1.6| ) and (1^) we have 



A(m)=ao + PST, 5(m) = oi log/i/m + 60 + PST , (1.9) 



and the FOMO approximation is given by 
da 



= aoaemtts + (ai log /i/m + 60) ctema^ ■ (1-10) 

FOMO 



dydp, 

Our final result will be given by 

FONLL = FO +(RS - FOMO) xG{m,p^). (1.11) 

The notation FONLL stands for fixed-order plus next-to-leading logs. Formula (|1.11| ) 
is our practical implementation of eq. ( |1.8D . 

The quantities appearing in the RHS of eq. ( |1 . 1 1| ) are available as Fortran com- 
puter codes. Actually, due to the non-physical splitting of photoproduction cross 
sections as given in eq. ( p. . 1|) , the pointlike and hadronic components are usually 
computed by different packages. It is therefore useful to write 

FO = FOpnt+FOhdr, FOM0 = FOM0p,t+FOM0hdr, RS = RSpnt+RShdr- (1.12) 

The package that computes FOhdr was originally developed in ref. P], and subse- 
quently modified in ref. |]T^, in order to implement FOMOhdr (whose analytical form 
was obtained in ref. [§]), and in order to use the appropriate renormalization and 
factorization scheme. FOpnt was computed in refs. 0] and 0, and FOMOpnt was 
obtained in the present work. The code for FOpnt (taken from ref. [|I|) has been 
extensively modified in the context of the present work in a manner analogous to 
what done for FOhdr in ref. [|^. The package that computes RShdr was presented 



in ref. |T^, suitably modified for heavy quark fragmentation in ref. W^, and used in 
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ref. [12|. Finally, the package relevant to RSpnt was written by the authors of ref. [14 
and adapted to heavy quark fragmentation in ref. 0. It has undergone some further 
modifications during the course of this work. It should be clear that FO, FOMO, 
and RS are all strictly interrelated. This fact provides us with a way not only to 
test that the various codes mentioned above are mutually consistent, but also that 
their implementation in our formalism, eq. ( |1.8|) , has been carried out correctly. In 
particular, the relation between FO and FOMO, and between FOMO and RS, will be 
the argument of sections ^ and ^ respectively. In the latter section, we shall show in 
particular that the FOMO and RS results differ only by terms of order aem«s- 

The paper is organized as follows. In sect. ^ we describe the procedure to adopt in 
order to translate the FO result from a scheme with rif — 1 light flavours to a scheme 
with rif light flavours. In sect. |^ we give a few details concerning the calculation of 
the massless limit of the FO calculation. In sect. § we check the matching between 
the FOMO and the RS calculation. Unlike in the hadroproduction case, subtleties 
arise here due to a different separation of the hadronic and pointlike contributions 
in the FO and RS approaches. We shall see that only in the full (i.e. pointlike plus 
hadronic) cross section we have complete matching. In sect, ^we examine the size of 
power-suppressed effects in order to understand at which value of m/p^ the massless 
approach gives a sensible approximation to the massive calculation. The function 
G{m,p^) will be chosen on the basis of the considerations given in this section. In 
sect, ^we describe our full result, for the case of charm production at HERA. Finally, 
in sect. |^ we give our conclusions. 



2 The change of scheme 



As shown in detail in ref. |T^, a change of scheme from the one with n\{ to the 
one with rif = riif + 1 light flavours brings about the following changes in and in 
the parton densities 



a("»)(^H,) = «f ^(/"K) - T^Tp log (^^) + o{at) 



371 



1 + 



■^og — 



(2.1) 



(2.2) 



All other parton densities are affected at higher orders in as. Eqs. ( p.l|) and (|2.2| ) 
are universal, that is, process-independent. It is easy to convince oneself that the 
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pointlike and hadronic components of the cross sections transform independently 
under this change of scheme. Thus, for what concerns the latter component, we can 



safely use the formulae of ref. |jT2|. In the pointlike component, the only effect at 
0{aemC(s) is generated by the Born- level 7(7 cross section. It is a matter of trivial 
algebra to conclude that, in order to go from the nif- flavour to the nf-fiavour scheme, 
the fixed-order cross section has to be modified by adding a term 

5a,, = -«,;lTplog4<^ (2.3) 

to the 7(7 cross section. For any reasonable range of scales, this correction is not large, 
and it vanishes for fip = /i^. 

In the following, we shall always refer to the FO and FOMO calculations performed 
in the nf-fiavour scheme. We shall thus always assume that as and the parton densities 
Fj refer to ai"'^ and 



3 Massless limit of the fixed-order calculation 

As in the case of the change of scheme, the massless limit of FO can clearly be 
performed independently for the pointlike and hadronic components. The latter has 
been considered in ref. |]T^; here, we only deal with the pointlike part. 

The massless limit of the fixed-order pointlike cross section formulae (in the sense 
of eq. (|1.10|) ) is obtained via algebraic methods from the results of ref. As pointed 



out in ref. [Q, the limiting procedure is non-trivial. In fact, the partonic cross sections 
at order aem«s contain distributions, such as delta functions or principal value sin- 
gularities. When taking the massless limit, new contributions to these distributions 
arise. We have computed this limit analytically by a computer-algebra procedure 
applied to the massive cross section formula, and checked its correctness in the fol- 
lowing way. We compute the heavy-flavour differential cross section at fixed p^, y, 
and centre-of-mass energy. We choose the renormalization and factorization scales 
equal to p^. Under these conditions, the mass dependence of the result is confined to 
the partonic cross sections. In the massless limit approximation, the only remnants 
of mass dependence are in logarithms of the mass in the 0{aem(^s) terms. Thus, if 
we plot the FOMOpnt cross section versus the logarithm of the mass, we get a straight 
line. On the other hand, if we plot the full FOpnt cross section versus the logarithm 
of the mass, it should approach the FOMOpnt result in the limit of small masses. We 
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Figure 1: Comparison of the FOpnt and FOMOpnt differential cross sections 
as a function of the logarithm of the mass, at = 20 GeV and y = 1. 



have performed this test, choosing the proton energy Ep = 820 GeV and the photon 
energy = 20 GeV. We adopt the CTEQ4M set [ffH] for the proton parton density 
functions. These parameters will be our reference choice from now on. The results 
are displayed in figs. |l] and ||. From the figures, it is quite apparent that the massless 
limit, as well as its implementation for the calculation of cross sections, was carried 
out correctly. There is also an important observation to make: the FOMOpnt cross 
section is larger than the massive calculation, i.e. power suppressed mass effects are 
negative, contrary to the case of hadroproduction. 

Notice also that the FOMOpnt approximation is quite accurate even at relatively 
large values of m/px- For example, from both figs. |I| and |] we notice that even for 
m/px ~ 1/2, the FOpnt cross section differs from FOMOpnt by at most 30%. This 
has to be contrasted with the case of the hadronic component (see ref. [0), where 
FOMOhdr is describing well FOhdr only for rather small values oim/p^. This mismatch 
between the behaviour of the pointlike and the hadronic component will clearly show 
up in phenomenologically relevant cases, as we shall later see in sect. 0. 
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Figure 2: As in figure for 73 (left) and (right) components alone. 



4 Matching 

We now examine the matching between the resummed approach and the FOMO 
calculation. In this case, there is a strict interplay between the pointlike and the 
hadronic components, that should be dealt with very carefully. 

There are ingredients in the resummed approach that are not explicitly present in 
the FOMO calculation. These are the fragmentation functions for final-state partons 
to go into the heavy quark, and the parton density for finding a heavy quark inside 
the hadron. The fragmentation function for any parton to go into a heavy quark has 
a power expansion in terms of the coupling constant evaluated at the scale yuQ, and 
of logarithms of /x/m: 

00 k 

D,(x,/i,m) = $^$^4'''H^) W-«s(/i), (4.1) 

k=0 1=0 

that can be obtained by solving the evolution equation for the fragmentation function 
at the NLL level, with the initial conditions of ref. [Q. Similarly, the parton density 
for finding the heavy flavour in a hadron can be expanded in the form 

00 k 

F,{x, f,,m) = J2Yl ^K/")) W - «s (/i) • (4.2) 

k=0 1=0 

With Filn) in the argument of the coefficients, we mean that the coefficients have a 
complicated functional dependence upon the parton densities evaluated at the scale 
^We take for simplicity /Ir /ip = Pt, and denote the common value with /i. 
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(a) (b) 
Figure 3: Photon splitting into a heavy-quark pair. 



/i. The existence of formal expansions of the form (^4.1| ) and (|4.2| ) can be easily 
proved, by writing the Altarelli-Parisi equations in integral form, and then solving 
them iteratively. A more detailed argument was given in Appendix A of ref. [0. 

Once eqs. (4J.) and ( |4.2D are formally substituted in the RS cross section formula, 
this formula itself becomes a power expansion of the form of eq. (|1.7|) , with the 
coefficients that depend (functionally) upon the structure functions for light partons, 
in the nf-flavours scheme, evaluated at the scale /i. The FOMO calculation has an 
expansion of the same form (truncated to order aemCtg) "with coefficients that are also 
dependent upon the same light-parton structure functions^ evaluated at the scale 
fi. Thus, because of the next-to-leading logarithmic accuracy of the resummed cross 



section, the terms up to the order 
calculation. 



in RS will match exactly with the FOMO 



In the photoproduction case, the matching can take place only in the full cross 
section, i.e. pointlike plus hadronic. In fact, let us consider the photoproduction 
subprocess in which a photon splits into a heavy-quark pair, and afterwards the 
heavy quark scatters with a parton coming from the hadron, as shown in fig. |^. The 
contributions depicted in the figure, in the kinematic region of a photon splitting 
coUinearly, is fully included in the pointlike contribution in the massive, fixed order 
calculation. In fact, because of the mass, no collinear subtractions are needed on the 
photon side. 

In the resummed approach, instead, these graphs are collinear divergent, and the 



^We observe that this property of the FOMO calculation is only valid in the modified scheme 
described in section ^ If we had used the standard scheme for the fixed-order calculation, the 
structure functions and the coupling as appearing there would be those with ni — 1 flavours. 
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singularities are subtracted at a scale /i. A corresponding hadronic contribution is 
present, where a heavy quark is found inside the incoming photon. Such a contribution 
is not present in the hadronic part of the fixed-order approach. 

Therefore, in order to check the matching, we should consider pointlike plus 
hadronic cross sections, and compute the difference between two relatively large num- 
bers (RS and FOMO), each of which is obtained as the sum of two numbers (RSpnt plus 
RShdr and FOMOpnt plus FOMOhdr)- In practice, this procedure requires an extremely 
careful treatment of numerics, which is not required anywhere else in our study. We 
thus adopted a slightly different, although equivalent, strategy, which is based upon 
the observation that, as far as the matching is concerned, it is only the contribution 
depicted in fig. ^ that is treated differently in RS and FOMO. Therefore, we can 
simply add to the resummed pointlike result the hadronic contribution with a heavy 
quark in a photon. Since the matching is checked up to the order Ue^a"^, it is enough 
to compute such contribution at this order. Calling s the s-channel invariant of the 
photon-hadron system and i/h iyi) the heavy quark (light recoiling parton) rapidity 
in the photon-hadron CM system, we immediately find 



dyhdp. 



.2 



j^X ^'..(-.) £| f^/-^'^^ dm (4.3) 



where PQ stands for photon to quark, and 

P,,{x) = '^N^el{x' + {l-xf). (4.4) 

d^hi /d^2 is the leading order (i.e. 0{al)) massless cross section (without the two- 
body phase space) for the partons h (the heavy quark) and I (a light quark or a gluon) 
to give an outgoing heavy quark h with transverse momentum and rapidity y in 
the parton centre-of mass (CM). The variables and y characterize completely the 
scattering kinematics in the partonic CM system. We have furthermore 

y = — 7, — . (4-5) 



X, = ^/M!^^pM±^^, (4.6) 



X = exp(-|/fe) +exp(-y;) ^^^^ 



and the partonic cross sections are 

rf$2 ^ 2s9 P 



'^S^^^-^^ - (4.a.)^ll!^, (4. 
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where the partonic Mandelstam variables s, i and u are easily obtained from p-p and 
y. For future reference, we shall call PQ (for Photon to Quark) the contribution of 
eq. ([4 .31 ). Such contribution is part of the direct component of the FOMO result. 

In order to check the matching of the resummed and fixed order calculation we 
proceeded in the following way: 

(a) we computed the FOMOpnt result in the Uf flavours scheme; 

(b) we computed the RS result without the intrinsic heavy quark component in the 
hadron parton densities, and with a heavy quark fragmentation function set 
equal to 6{1 — z) for the Dg component, and all other components set to zero; 

(c) we computed the LO contribution to the RS result, with the heavy quark com- 
ponent in the hadron parton densities set equal to 

/TT _l ^ Z 



and all other components set to zero; 

(d) we computed the LO contribution to the RS result, with the heavy quark density 
set to zero, and the heavy quark fragmentation function set to its 0{as) value 



27T 



^og4-21og(l-^)-l 



(4.11) 



Di{z,fi) 0, for i^g,h. (4.13) 



D,(z,/.) ^f£^(.^ + (l-,)2)logi^ (4.12) 



We then verified that the sum of items (b), (c) and (d), plus the PQ contribution 
(eq. ( |4.3| )) is exactly equal to item (a). Observe that with this procedure we were 



able to isolate the terms of order aemtts and aemOs the RS result, and thus, in 
order to check the matching, we did not need to go to the weak coupling limit, as was 



done in ref. [T^. The resummed hadronic component, suitably augmented of the PQ 
contribution, must then match the fixed order hadronic component in the massless 
limit. This follows easily from the work of ref. ]T^|. 
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Figure 4: FOMOpnt compared to the sum of RSpnt and PQ. The parameters 
are as in fig. |l]. The charm quark mass was set equal to 1.5 GeV. 



The procedure we carried out, besides convincing us of the correctness of our 
theoretical approach, has also served as a test of consistency between the computer 
programs of ref. and the fixed order calculation of ref. 

We have thus demonstrated that the quantity RSpnt + PQ — FOMOpnt is of order 
aem«s • However, this does not guarantee that it is also small in practice. In ref. [|12 



it was shownQ that RS - FOMO in the case of hadronic collisions is of order , but it 
is numerically non-negligible even for bottom production. We shall now show that, in 
the present case, the difference RSpnt + PQ — FOMOpnt is numerically small, at least in 
the region where Oslog Pt^/itl is also small (that is to say when is not large). This 
is in fact what we see in fig. |[ We also see that, as previously discussed, FOMOpnt 
is not a good approximation of RSpnt, due to the lack of the photon-splitting term 
■^See in particular fig. 8 of ref. [O 
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Figure 5: FOpnt and FOMOpnt at Born and full O(aemas) level, plotted 
as a function of the mass and at fixed transverse momentum (left figure) or 
fixed transverse mass (right figure). 



PQ. We further note that the matching deteriorates at higher values of due to the 
increasing importance of the resummation performed in RSpnt but absent in FOMOpnt- 

5 Power effects in the RS and FOMO calculations 

When comparing and matching the FO, FOMO and RS approaches, there is much 
arbitrariness in the way mass effects are treated. For example, we may decide to 
compare transverse-mass distributions instead of transverse momenta. These are 



equal for the FOMO and RS calculations, but differ in the FO case. In ref. |I2| the 
effect of this replacement was studied for hadron-hadron cross sections. Here we thus 
study only the pointlike component FOpnt and FOMOpnt- 

In fig. H we plot FOpnt and FOMOpnt as a function of the mass, keeping either 
or fixed. It is clear that, in the plots at fixed transverse momentum, the massive 
calculation is more suppressed near the threshold (i.e., as m approaches p^)- In this 
region higher Xp values are probed in the proton structure function with respect to 
the massless calculation since, besides the transverse momentum, also the mass has to 
be produced. This is clearly a spurious effect, and we shall always prefer to perform 
the matching at fixed transverse masses. The right plot of fig. ^ is an a posteriori 
justification of this procedure, since it is quite clear that in this case the massive and 
massless limit results are much closer to each other. It is also worth noting that in the 
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Figure 6: As in the right plot of fig. ^, for y = —0.6, 2 and 3. 



right plot of fig. ^ power suppressed mass effects are positive, in contrast to what we 
observed earlier (see fig. |l|). In figure |^ we show analogous plots at fixed transverse 
mass for y = —0.6, 2 and 3. In all cases the massless limit and full massive results 
are in good agreement, although a worsening is seen for the case of y = —0.6. This 
rapidity value is quite close to the phase space boundary, where we do in fact expect 
some anomalous behaviour due to Sudakov logarithms. 

We shall now proceed as follows. For a given transverse momentum the FO 
cross section is evaluated and combined, using eq. (|1.11| ), with the FOMO and RS 
results evaluated at the corresponding = a/Pt + value. In this way the three 
calculations are performed at the same m^. Moreover, a central choice for the renor- 
malization and factorization scales will be /Ir = /Xp = A/p^T-pm?, so that they coincide 
in the three calculations. 
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In the hadroproduction case |]T2|, a suppression factor G{m,prj,) was introduced to 
multiply the RS — FOMO term in eq. ( |1 . 1 1 ) , with G{m,pT,) approaching one at large 
Pt- The following form was chosen: 



Pt + c^m"^ 



(5.1) 



Adopting the same form, our final formula becomes 



FONLL = FO + 



Pt + C 



■[RS - FOMO] . 



(5.2) 



In ref. ||T2[ c = 5 was chosen, in order to suppress meaningless large corrections coming 
from the massless approach at low momenta. A detailed discussion of the role of c in 
the present case will be given in the next section. 



6 Phenomenological results 



We now present some benchmark results of our calculation, based upon eq. ( |5.2| ). 
In order to perform the calculation, we have suitably modified the codes relevant to 
FOpnt, FOMOpnt, and RSpnt, in order to obtain a consistent implementation. The 
hadronic quantities, on the other hand, have been produced with the program of 



ref. [0, here modified in order to allow the use of the photon parton densities. 

When applying our matched formalism to phenomenology, we must first make 
sure that the parton densities we use correctly incorporate the logarithms of m/fi 
that we are trying to resum. This must be the case if the structure functions have 
been evolved correctly, with the heavy flavour evolution turned on when /i = m, as 



appropriate in the MS scheme |]T6| . In the structure function fits available today this 
is not always the case. In part this is due to the fact that the parton densities are 
often given as an interpolating grid, and the region near the heavy flavour threshold 
may not be represented accurately. Furthermore some parton density sets do not 
implement the heavy flavour thresholds according to ref. |T^. While for hadron 
structure functions several sets with a correct charm density are available, the choice 
among photon structure functions is quite limited. We have chosen the AFG set [0 , 
which is in the MS scheme, and claims a correct implementation of the charm density. 
In order to test this fact, we have plotted in flg. |^ the AFG charm parton density in 



16 




the photon together with a charm density computed with the following formula: 



271 



2tt 



dz 



)Fg{x/z) — 



1 fj'l 

log^ 



(6.1) 

(where Cc = 2/3 is the electric charge of the charm quark) that should hold for 
factorization scales not too far from the heavy quark mass. In order to match the 
AFG parameters, we have chosen m = 1.41 GeV in this plot. We expect that the 
slope of the AFG charm density and our approximate formula should agree for small 
scales. Unfortunately we do not observe this, especially at small values of x. We 
have checked, however, that these differences would only be marginally important if 
our approximate formula, eq. (|6.1|), were used to calculate physical cross sections, the 
results obtained with the original AFG charm densities being well reproduced. 
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Figure 8: Contours in tlie y-pT plane of constant minimum value probed 
in the hadronic component (on the left), and of constant ratio of the direct 
component over the total cross section (on the right) for the NLO result. 
The short-dashed line represents the phase space boundary. 



In what follows, we consider our reference case of photon-proton collisions, with 
= 20 GeV and Ep = 820 GeV. In the left panel of fig. ^ we plot contour lines 
representing the minimum value of allowed in the hadronic component of the cross 
section, in the y-pT plane. Commonly available photon structure function fits are 
undefined below values of x-y of the order of 10~^. The minimum value of x^ = 1.5 10^^ 
relevant to the AFG set is explicitly shown in the figure. For rapidities currently 
probed at HERA, smaller x^ values are never reached. In the right panel of fig. ^ 
we show, for future reference, the ratio R of the pointlike component over the total 
cross section, both computed at the NLO. Where the hadronic component prevails 
(i.e. R becomes small), we expect to find the same problems found in ref. as far 
as the matching between RS and FOMO is concerned. We point out that the results 
displayed in the right panel of fig. || have a scale dependence of 0{aemC(^), and they 
are thus accurate to LO only; however, the plot gives a clear idea on the dominance 
of pointlike or hadronic components in the physical cross sections. 

In our phenomenological study we shall use m = 1.5 GeV, the AFG set for the 
photon and CTEQ4M for the proton parton densities. In fig. ^ we show the rapidity 
distribution for charm quarks at = 2 GeV. In the left plot of this figure, both the 
fixed order prediction and our FONLL results are shown. The photon has negative 
rapidity in our convention. The phase space limit for the rapidity ai p^ = 2 GeV is 
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Figure 9: Rapidity distribution of charm quarks at = 2 GeV, with no 
smearing factor for the term RS — FOMO (i.e., c = 0). In the right figures 
the various contributions to the FONLL cross section are shown. 



—2.77 < y < 6.49. Thus, the left limit of the plots is near the negative rapidity limit. 
No suppression factor was applied to the difference RS — FOMO. At these low values 
of Pt, we would expect that the fixed order and the FONLL result should agree. In 
fact we observe good agreement at negative rapidity, while at moderate and large 
rapidity the FONLL result is much larger than the fixed order one. 

A detailed analysis of the various contributions to the FONLL cross section is 
shown in the right plot of the figure. We show the RS result, and the mass corrections 
FO-FOMO, for the pointlike and hadronic component separately. For comparison, also 
the fixed order pointlike and hadronic component are shown. 

A few comments are in order. We begin by noticing that the difference FOpnt — 
FOMOpnt (the dashed line in the left lower corner of the plot) is quite small compared 
to FOpnt. On the contrary, FOhdr — FOMOhdr is even larger than the FOhdr result, 
due to a negative FOMOhdr contribution. Comparing the pointlike and hadronic fixed 
order results, we see that the first one prevails up to rapidities of 0.4. A similar plot. 



for Pt = 3 GeV, is shown in fig. |10|. We see similar features to the case of fig. |[ 

We conclude that the large difference between our FONLL and FO results is only 
due to the hadronic component of the cross section. 



This is consistent with what already observed in ref. [12|. In that paper these 



differences were ascribed to a large RS — FOMO term which, due to spurious higher 
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order terms present in RS, did not vanish quickly enough in the small Pt limit. 
Consequently, the small-pT suppression factor G{m,pT) = p\/ {j)\ + c^rn?), with c = 5, 
was applied to the difference RS — FOMO in order to get rid of this problem. 

In this work we can now make a conjecture about a possible relation between 
the higher order behaviour of the "massless" RS calculation and the size of the mass 
terms: in the hadronic case mass terms are large and so are the higher orders in RS. 
In the pointlike component, instead, the mass terms are relatively small, and there 
appears to be a good cancellation between RS and FOMO. One can therefore make 
the reasonable assumption that the less important the power suppressed mass terms 
are, the better behaved a resummed "massless" approximation will be. 

Having said so, we observe that for the pointlike component alone no suppression 
factor is actually needed, consistently with what we observed in section ^. 

At this point, one is tempted to apply the suppression factor G{m,prj,) only to the 
difference RShdr — FOMOhdr- However, this would be incorrect. As we have already 
discussed in section H, the terms of order aem«s and aem«s in RShdr do not match 
FOMOhdr- The difference RShdr — FOMOhdr is of order aemC^s- Thus, using a different c 
value for the pointlike and hadronic component would lead to the introduction of mass 
suppressed terms of order aem«s5 ^^at would spoil the accuracy of our calculation. It 
is instead sensible to use different c values in the following expression 

2 

FONLL = FOp„t + FOhdr + (RSp,t-FOM0p„t + PQ)^-^^ ^ + 

Pt + Cpnt^ 
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Figure 11: Effect of the smearing functions in the FONLL result (eq. (|6.2|) ) 
for the distribution of charm quarks at y = 0. In the right figure, the 
same results are multiplied by m^- 



(RShdr - FOMOhdr - PQ)^^i ^, (6.2) 

since now the expressions in parentheses are indeed of the appropriate order. From 
formula ( |6.2| ) with Cpnt = Chdr = c we recover the standard expression eq. ( |5.2| ). In 
fig. |ll] we show the effect of varying independently Cpnt and Chdr in the FONLL result 
of eq. ( |6.2| ). It is quite clear that the inclusion of a smearing function for the pointlike 
component has little effect (the difference between the two solid lines). The smearing 
function for the hadronic component has instead a very large effect, suppressing the 
cross section at small p^. In its absence (dashed curve), one gets a cross section that is 
roughly twice as large as the fixed order one. At large Pt the sensitivity to a smearing 
function decreases rapidly. The FONLL approach gives smaller cross sections than 
the fixed order approach in this region. This is what one expects from resummation 
effects, since the emission of collinear gluons, in general, has the effect of softening the 
spectrum. Similar plots are also shown in fig. |l2|for y = —1. Also here we see that 
the effect of the smearing function on the pointlike component of the cross section is 
small, and we see a relatively large effect of the smearing in the hadronic component. 
This effect is not quite as large as in the case of y = 0, because at negative rapidity 
the hadronic component is less important. 

We shall now show a few benchmark results for photoproduction of heavy flavour 
at HERA. For simplicity, we shall always keep c = Cpnt = Chdr = 5. The aim of these 
results is to assess the difference between the fixed order results, the resummed results, 
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Figure 12: As in fig. |Tl|, for y 
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Figure 13: Comparison between the FONLL, FO and RS results. 



and our matched formahsm. We assume an incident photon energy of 16.5 GeV 
(corresponding to an electron energy of 27.5 GeV and to a photon energy fraction of 
0.6), and a proton energy of 820 GeV. In fig. |1^ we show a comparison between the 
FONLL, FO and RS results. We see that the RS and FO results are remarkably close 
for small Pt, a fact that we have understood to be true for the pointlike component, 
and purely accidental for the hadronic component. At moderate Pt also the FONLL 
result is very close to the FO result. While in the case of the pointlike component 
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this is true in all cases, for the hadronic component this only happens if one chooses a 
relatively large value for c. At very small the resummed calculation becomes totally 
unreliable. At large p^, the FONLL and RS results become quite similar (since, by 
construction, the G{m,prj,) suppression factor tends to one, and FO tends to FOMO, 
canceling it - see eqs. ( |1 . 1 1|) and (|5.2|) ), and remain smaller than the FO result. This 
is a general consequence of the effect of multi-gluon emission resummed in the RS 



and FONLL approaches. The pattern displayed in fig. 13 seems to be quite universal. 



Similar plots for y = 1 and y = —1 are shown in fig. Furthermore in fig. [Ig we 
also show the case of photon energy fractions of 0.4 and 0.2 (that is -E^ = 11 and 5.5 
GeV respectively). 



7 Conclusions 



In this work we have implemented a technique to compute the transverse momen- 
tum spectrum in heavy flavour photoproduction, which is accurate to the full NLO 
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level at moderate values, and to the NLL level at large p^- 

This is achieved by properly merging a "massless" resummed approach, valid in 
the large-px region, with a full massive fixed order calculation, reliable in the small-pT 
one. 

We observe that, for the pointlike component of the cross section (in the sense 
discussed in section ^), the massless limit is a good approximation to the full cross 
section, provided one uses ttIt rather than in the former. On the contrary, this 
is not the case for the hadronic component of the cross section, as already observed 
in ref. []T2[ in the hadropro duct ion case: power suppressed mass terms happen to 



be important in the moderate region. This suggests that a proper merging like 
the one studied in this paper is necessary and superior to just employing a massless 
approach. 

The results obtained with our procedure are in good agreement with the fixed- 
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order calculations at moderate transverse momenta, and with the so-called massless 
resummation (or fragmentation function) approach at very large momenta, effectively 
interpolating, in a theoretically sound manner, between the two approaches. 

In general, we find that inclusion of resummation effects brings about a softening 
of the Pt spectrum at large transverse momenta. 
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